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Abstract 

Fokker-Planck equations have been applied in the past to field theory top- 
ics such as the stochastic quantization and the stabilization of bottomless ac- 
tion theories. In this paper we give another application of the FP-techniques 
in a way appropriate to the study of the ground state, the excited states 
and the critical behaviour of quantum lattice Hamiltonians. With this pur- 
pose, we start by considering a discrete or lattice version of the standard FP- 
Hamiltonian. The well known exponential ansatz for the ground state wave 
functional becomes in our case an exponential "cluster" expansion. With a 
convenient choice for this latter, we are able to construct FP-Hamiltonians 
which to a large extent reproduce critical properties of "realistic" quantum 
lattice Hamiltonians, as the one of the Ising model in a transverse field (ITF). 
In one dimension, this statement is made manifest by proving that the FP- 
Hamiltonian we built up belongs to the same universality class as the standard 
ITF model or, equivalently, the 2D-classical Ising model. To this respect, some 
considerations concerning higher dimensional ITF models are outlined. 
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1 Introduction 



A fundamental problem in Statistical Mechanics and Solid State Physics is the search 
of the ground state and excitations of many body Hamiltonians. Exact solutions 
have been obtained for integrable one dimensional Hamiltonians PJ,P|], [[§• However, 
the majority of the models are not integrable, specially in dimensions greater than 
one, and in these cases one has to resort to approximation techniques. We have 
squematically, 



H (Hamiltonian) 



Perturbation theory 
Variational methods 
Renormalization group 
Numerical approaches 



ipo) (ground state) (1.1) 



The path from H to | -00 ) is very cumbersome and can only be fully accomplished 
in a few cases. Instead, one could try to reverse the arrows in (|1 . 1| ) , by making a 
clever choice of the state | ipo) first, and asking latter for the Hamiltonian H for 
which the proposed ansatz becomes the exact ground state (this second strategy 
can be paraphrased in the words of Q as "finding questions to some interesting 
answers". Some remarkable cases where this "inverse method" has been applied 
would include the exponential Jastrow wave functions, used in the study of liquid 
4 He and nuclear matter ||, the valence-bond solid states (VBS) 0, the Laughlin 
wave function M, etc. The path from | ?/>o) to H could be summarized as follows: 








[1.2) 



where a denotes either sites in a lattice, or particles or some other quantum numbers. 
In words, one has first to find a set of operators R a which annihilate the state | V'o); 
afterwards one builts up a positive semidefinite Hamiltonian H according to the 
prescription ( |1.2|) . If the state | ipo) of ( |1.2|) is a valence-bond state, then the R a 
are projector operators (i.e. R^, = R a , R^ = R a ) and we can write H = J2 a Ra- 
However, for Jastrow type wave functions, the R a are not projectors, and these will 
be the kind of operators which we shall be concerned with. 

The aim of this paper is to combine the direct and inverse methods ( |1.1| )-( |IT2"| ) in 
the study of quantum lattice Hamiltonians. To be more concrete, the problem we 
want to adress is the construction of the ground state | ipo) and the excited states of 
a quantum Hamiltonian H defined on a lattice. The strategy could be summarized 
in the following diagram: 



Perturbation Fokker — Planck 

theory method 
H — > ip (Exponential ansatz) — > H (1-3) 



The first step consists in the construction of a state of exponential type, which 
will be our candidate for the ground state of H. To make this ansatz, we shall 
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follow the perturbative-variational procedure of determination of the ground state 
| ^q(A)) of a lattice Hamiltonian, in the form H(X) = Hq + XHi, established in ||; 
here Hq is a Hamiltonian whose non degenerated ground state | 0) is known, and A 
is a coupling constant. 

Two alternative ways for writting the state | ipo(\)) of H(X) were proposed in 
reference ||, to wit: 

| Vo(A)> = exp(£ xnu n) I 0) (1-4) 

n>l 

or 



|^o(A))=exp(^a / (A)V7) | 0) (1.5) 
i 

where both U n and Vi are a certain set of interrelated local operators (see below), 
whereas the aj are certain functions of A. 

Expansion ( |1.4| ) is adequate to make contact with the standard perturbation 
theory (PT), whereas expansion ( |1.5|) is a kind of cluster expansion. Let us describe 
in some more detail the main ideas of M, since this paper is, after all, an extension 
of that work. 

When one expands the exponential of ( |1.4| ) and compares the result with the 
usual PT, one obtains a set of equations for the operators U n . Thus, we can say 
that PT to order v fixes , up to some freedom, the familly of operators {U n }\< n < u . 
Working through particular examples, as the Ising model in a transverse field, one 
observes that each operator U n can in fact be written as a sum of "irreducible" 
operators Vf. 

U n = Y,Pn,lVj (1.6) 

where in || the label / denoted a cluster formed by a finite number of l's distributed 
on certain sites of the lattice, the rest of sites being occupied by zeros. The operator 
Vi obeys a fusion algebra given by VjVj = Vj + j, in such a way that the cluster I + J 
is obtained through the usual Z 2 -sum rules: + = 1 + 1 = ,0 + 1 = 1 + = 
1 (this rule has its origin in the Z2-symmetry of the Ising model). The usefulness 
of eq.( |1.6|) lies in the fact that a given cluster operator Vj may appear at different 
-and possible at all- orders in PT. This would imply that for a fixed value of /, p n i 



might have an infinite number of non zero values. Introducing ( |1.6|) in (|1.4j ) we can 
relate the perturbative and cluster expansions with: 



«/(A) = E Pn.A n (1-7) 
n>ui 



where the number 



vj = minimun{n/p ni 7- ^ 0} (1-8) 

is called the level of the operator Vi\ somehow vj measures the size of the cluster 
0. Definition ( |1.7| ) implies that ai (A) has a Taylor expansion starting at X UI and 
containing non vanishing terms of higher powers whenever Vi appears in U n . 
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If we think of matrix with two entries n and / -or rather i/j-, we see by- 

its definition that p n> i is a lower triangular matrix of infinite dimension (see fig. la). 
Thus PT to order v gives all the entries of p n j to a "depth" v, i.e {p n ,i}i<n<u (see 
fig. lb). 

However, eg- flOp provides another approach to the exact ground state. Indeed, 
the knowledge of a single function «r(A) for a given operator Vi involves all orders 
in PT. This result suggests a non-perturbative approach to | ip ) based on the 
truncation of ( |1.5| ) to just those operators Vj with an order less or equal to a given 
order v (see fig.lc): 

|^>=exp( ]T o>V 7 )|0> (1.9) 

J, UJ<U 

In || the functions ocf^ were fixed by using the state | i) {u) ) trial variational 
state for the Hamiltonian H(X) = Hq + \H\. In this case the variational ansatz 
agrees to order v with the perturbative ground state and gives non perturbative 
predictions for large values of A, where PT cannot be trusted. 

In this paper we propose an alternative way for computing the parameters otf\ 
The principle of our method is a generalization of the well known Fokker-Planck's 
procedure, by means of which one associates a Hamiltonian to a wave function having 
an exponential form (see appendix A, for an application to Quantum Mechanics). 
The rule is simply to act with H Q on the proposal ( |1.9[ ) until one arrives to a Fokker- 
Planck Hamiltonian H FP '^{ocf^) which depends, of course, on the whole set of 
parameters ocj entering in the definition of | ip^). The condition we use to fix 
these parameters as functions of A is that H FP '^ should only differ from H in 
operators of a level greater than v, i.e: 

H-H FP ^= Cf\\)V! (1.10) 

I, vi>y 

Correspondingly, the FP-energy E FP ^ V \\) although incorporating non pertur- 
bative effects, should agree with the exact ground state energy up to order \ u+1 . 
Compairing E FP,U (X) (for v = 1 and 2) with the variational energy E var,u (X) ob- 
tained in || in the study of the Ising model in a transverse field, one observes that 
the latter gives better values than the former, as to the energy is concerned. It 
seems natural then to wonder about the advantages of the FP-approach versus the 
variational one. Well, the problem is that within the framework of the variational 
approach is very difficult to go to higher values of v since, for example, the norm of 
the variational state ( |1.9j ) turns out to be the partition function of a statistical me- 
chanical model whose couplings are dictated by the operators Vj. On the contrary, 
in the FP procedure the norm of the state ( |1.9| ) is not needed to fix the parameters 
a j in terms of A; rather, what is required is the algebra satisfied by the cluster 
operators Vi, and this algebra is relatively simple in the example considered in this 
paper. Besides, within the FP-method we can always match the perturbative re- 
sults to any order, at the expense of including a sufficient number of operators in 
the ansatz ( |1.9|) . It is clear, however, that except for the lowest orders in PT the 
complexity of the ansatz will obscure its usefulness. On the other hand, one of the 
aims in constructing more general ansatzs for ground states is to see whether or not 
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they exhibit critical behaviour. In this sense, what actually matters is to determine 
the universality class which this Hamiltonians belongs to, rather than their appar- 
ent simplicity. Bearing in mind these ideas, we can try to reinterpret eq.( |1.3|) as a 
map from a given Hamiltonian if to a FP-Hamiltonian H FP belonging to the same 
universality class as H; in this case H and H FP would share the same critical prop- 
erties. To achieve this map one should consider ansatzs involving cluster operators 
Vi with all the possible sizes (see fig. Id). We shall show in this paper an example 
where such a construction can be accomplished. 



The organization of the paper is as follows. In section 2 we review the application 
of the FP machinery to field theory in the continuum and reconsider it in such a 
way that it can be applied to lattice Hamiltonians. In section 3 we benefit from 
this generalization to look into the Ising model in a transverse field (ITF) in D = 
1 and D = 2 dimensions. In section 4 we generalize the exponential ansatz of 
section 3 in order to include operators at all length scales, and we show in the 
one dimensional case the existence of a critical behaviour for a particular value of 
the coupling constant. In section 5 we study the problem of the excited states of 
FP-Hamiltonians using a "change of variables" which brings the Hamiltonian to a 
simpler form H RFP , which we call "reduced" FP-Hamiltonian. Finally, a general 
discussion with prospects for the future is proposed in section 6. 

We have also included four appendices which illustrate in detail some aspects 
mentioned in the course of the paper. In particular we collect in appendix D the 
results of sections 3 and 4 extended to the case of ITF models in dimensions higher 
than one. 



2 Fokker-Planck lattice hamiltonians: 
General construction 

The Fokker Planck Hamiltonians arise naturally in the study of stochastics processes 
in Statistical Mechanics or Quantum Field Theories. The idea is to consider a 
particle or a field evolving in time under the combined action of a force and a 



gaussian noise fllU| . The evolution equation for such a field 4>(x,t) is given by the 

+ Tj(x,t) (2.1) 



Langevin equation: 

d(j)(x,t) 5S[(j)] 



dt 6, 

where S[<j)} is the potential and r](x, t) is the noise. Here x G R D is the space location 
of the field. The corresponding probability P[(j), t] for finding a given configuration 
4>(x) at the time t is governed by the Fokker-Planck (FP) equation. A suggestive 
way of writing this equation comes up if we make the change of variables P[4>,t] = 
e ~ -0[0,t]; with this change, the FP equation for ip[(f),t] becomes a Schrodinger 
like equation in imaginary time, namely: 

^ = -*™VbM (2-2) 

where 
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is called the Fokker-Planck Hamiltonian. To this point, several comments are in 
order: 

(i) H FP is a positive semidefmite operator. This follows inmmediately if we write 

the Hamiltonian H FP as : 

H FP = i J d D xR\x)R{x) (2.4) 

where 

^ =i m +i 2m (2 ' 5) 

(ii) The functional ipo[<j>] = e~~^ is annihilated by the family of commuting oper- 

ators R(x) (i.e. R(x)R(y) = R(y)R(x) ) 

R(x)ip = 0, Vx (2.6) 

Hence, if ipo is a normalized state, then it becomes a ground state of the 
FP-Hamiltonian (O). 



(iii) The norm of the ground state ipo [0] can be interpreted as the partition function 
of a statistical mechanical model or euclidean quantum field theory in D- 
dimensions, with Boltzmann weights e" 5 ^' 



(Vo I -0o> = J D<P e~ s M (2.7) 



These simple properties have been used in a variety of applications, some of 
which are worth mentioning. To begin with, let us quote the stochastic quantiza- 



tion method of Parisi and Wu [11]. Introducing a ficticious time (the fith-time), 
this procedure enables one to compute quantum mechanical expectation values of 
observables Q[<f>] in terms of time averages of fields evolving under the Langevin 
equation (|2.1|): 

(Q> = i / D<p Q[<j>] e- s ^ = lim I F ' Q[<j>(x,t)]dt (2.8) 

Z J T->oo 1 Jo 

This formulation has been applied to quantum field theories, including gauge 
theories. 

Another fruitful application carried out in reference [Q concerns to the stabi- 
lization of bottomless field theories, i.e. theories with a potential unbounded from 
below and, therefore, with divergent partition functions. Here one uses the un- 
bounded action S[<j>] and eq. ( ]2.3| ) to construct a FP-Hamiltonian, which thanks to 
property (i) is well defined. Since H FP is bounded from below, the corresponding 
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ground state i/jq 07 " will be normalizable; that gives rise to a bounded effective action 

S e ff[<t>] 

S e ff[<p\ through the definition tpQ OT = e 2 . This stabilization procedure has also 
been applied to a number of theories including gravity and 2D-matrix models 

g, 



Finally, we want to mention the connection between FP-Hamiltonians and super- 
symmetry Let us suppose that we enlarge the field theory presented earlier with a 
pair of anticommuting fields 9(x) and 9(y) satisfying the anticommutation relations: 

{9(x),9(y)} = {9(x),9(y)} = (2.9) 

{e(x),6(y)} = S(x-y) (2.10) 
Now we define two super symmetric charges Q and Q 

Q = J dx R{x)9{x) (2.11) 

Q = Q ] = J dx R\x)6(x) (2.12) 

and satisfying 

Q 2 = Q 2 = (2.13) 

Then it is easy to see that the eq. ( p.6|) characterizing the state ipo can be trans- 
lated into the equations 

Wo = Wo = (2.14) 

In terms of Q and Q one can define a supersymmetric Hamiltonian Hss whose 
restriction to the bosonic sector (i.e 9 = 0) is nothing else but the FP-Hamiltonian 
73]), namely: 

Hss = \{QQ + QQ) 
= H FP + l -Jd D x -^[9(x),9(x)\ (2.15) 

S[4>] 

In the SUSY-terminology we can say that the exponential state tpo = e 2 is a 



supersymmetric ground state f!6 



The connection between the FP-method and SUSY is very suggestive but we 
shall not pursued it here, except for a discrete version of eqs.(|2.9| ) - (|2.15|) . 



As for the use of FP-ideas we make in this paper, let us say that they can 
be simply illustrated in the continuum as follows. Suppose we are given a field 
theoretical hamiltonian: 

H = jd D x[- l -^ 2 +Um (2-16) 
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where U{4>) is a funtion of <f)(x), V0(x), etc. Then if we look for the ground state of 
H in the form ipo[4>] = e~~^ , it turns out that the potential U(<f>) and the ground 
state energy E are related to S(4>) by the eq. 



Eq. ( p,17| ) is a functional version of the Ricatti equation | I7| . Consequently, the 
eigenvalue problem for the ground state of H is reduced, assuming the exponential 
ansatz V'oM = e ~ , to the solution of the functional Ricatti equation ( |2.17| ). 

In the context of Quantum Mechanics the analogue of eq. (|2.17 ) can be taken as 



the starting point of a perturbative- variational method []17| very close in spirit to the 
one developed in [[J. There are situations where the field theoretical Hamiltonian 
fl2.16|) can be diagonalized exactly, for example in terms of creation and annihilation 



operators as in pure QED |Tj|. The corresponding wave functional if>o[<f>] is then 



essentially a gaussian. When interactions are present, matters are not so simple. 
In any case, before looking for the ground state of quantum field Hamiltonians 
defined in the continuum, one should first regularize the model so as to control the 
infinities and this can be done, for example, by using a lattice regularization. This is, 
therefore, another motivation to consider the FP-method in the context of quantum 
lattice Hamiltonians. 



In the remaining of this section we shall present a generalization of the FP tech- 
niques introduced above. From an abstract point of view, the basic elements required 
in a FP-construction are: a) an algebra of operators, b) a vacuum state | 0), c) a 
free Hamiltonian H and d) an action S. In the continuum FP-construction, these 
elements are given respectively by the operators and <f>(x), the free Hamiltonian 

H = — | J d D x m&fp wnose ground state is the constant wave function, and the 
action S[4>] which is a functional of the field 4>(x). 

The construction can be implemented as follows: 



a) Algebra of operators. We consider a family of unitary operators {X a } ae c , 
{ZaIacC! satisfying the relations: 

X a X b = X b X a 

ZaZb = ZbZa (2-18) 
X q Za = uj a AZ A X a 

where the ui a A are phases (| u a A \— 1). The labels a and A run, respectively, over 
discrete sets C$ and Ci, whose nature will be specific for each model. In practical 
examples they will denote sites, links, plaquettes or internal quantum numbers. 

We shall demand that if Co and C\ contain an operator, say X a or Za, then they 
must also contain their hermitean conjugated, i.e: 

X\ = X a , a, a e C 
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Z A = Z A , A,AeCi (2.19) 

The Hilbert space where the operators X a and Z A are acting will be denoted by 
Ji. The algebra ( 2.18j) is actually the algebra satisfied by the order-disorder variables 



of, for example, the Ising model, the Potts model or the Z N gauge models. 

b) Vacuum vector. We shall suppose that in the Hilbert space 7i there is a 
vacuum vector | 0) satisfying: 

X a | 0) =| 0) , Va G C (2.20) 

The existence of a common eigenvector for all the X a -operators is in principle 
possible, because all these operators commute among themselves. Nevertheless, this 
vector needs not to be unique, i.e. there may exist other vectors satisfying Q2.20D ; 
this will be the case of models with a broken phase. 

c) Free Hamiltonian Hq. Let us define the following positive semi definite 
Hamiltonian: 

#o = ^E (*t - 1) (X a - 1) (2.21) 

which, using the unitarity properties of the X a operators and eq. (|2.19|) transforms 
in 

#o = £ (1 - X a ) (2.22) 

oeC 

With this definition, it is clear that the vacuum state | 0) is a ground state of H 

H | 0) = (2.23) 

d) Exponential ansatz. Suppose we are given a generic operator S(Zi, Z2, ■ ■ ■) = 
5(^4], which we shall call the action operator, and constructed out of the family of 
commuting Z-operators. From S[Za] we introduce the state: 

I ^ ) = e ~¥\- z A\ 1 0) (2.24) 
It is easy to check that the action of an operator X a on | 0) is given by: 

X a I Vo) = e-^^Xa I 0) 
= e -\s[u aA z A \ e \s\z A ] J ^ = e A a (s) 1 ^ ^.25) 



where we have defined 



MS) = \{S[Z A ] - S[u aA Z A }) (2.26) 



For latter purposes we shall write eq. ( |2.25| ) as 
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Ra | ^O) = 

R a = e- Aa(5) X a - 1 (2.27) 
Incidentally, it is noteworthy the commutation relation for the R operators 

[R a ,R b ] = (2.28) 

At this moment we may wonder whether there is a Hamiltonian for which the 
state | ipo) becomes an exact ground state. There are several answers to this question 
but in this paper we propose the following: 

Fokker-Planck Hamiltonian. Let Q a be a hermitian (fit — an d positive 
semidefmite operator (Q a > 0); then we shall call a FP-Hamiltonian the following 
hermitian positive semidefmite operator: 

H FP = \ E RfaaRa (2.29) 

Obviously the state ( |2.24| ) is a ground state for this Hamiltonian 

H FP | Vo) = (2.30) 

Unfortunately, the choice Q a — 1 although natural, gives rise to Hamiltonians 
where the operators X a and are mixed together. Yet, a peculiarity of this choice is 
that the corresponding FP-Hamiltonian presents a hidden supersymmetry. Follow- 
ing the same steps as in the continuum, this suppersymmetry can be made manifest 
if we introduce a discrete set of anticommuting operators 9 a and Ob satisfying: 

{0 a ,0 b } = {0 a J b } = (2.31) 

{0 a ,0 b } = 5 ab (2.32) 
The supersymmetric generators analogues to eqs.(|2.11|)-(|2.12j) read now 



Q= E #A ( 2 - 33 ) 

aeCo 

Q = E R % (2.34) 

aeCo 



and verify Q 2 = Q 2 = 0. 

The SUSY FP-Hamiltonian is then 



Hii = \{QQ + QQ) 
l -Y.RlRA + \Y.[ R a,R\]e a o a (2.35) 
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Hence in the bosonic sector Hgg reduces to eq.( |2.29| ) with Q a = 1. The problem 



with the Hamiltonian ( |2.35|) is that it mixes the X and Z operators in a complicated 
form, namely: 

<=i = E cosh(A a (S)) (e A ^ - X a ) (2.36) 

aeCo 

For this reason, and in spite of its intrinsic theoretical interest, we shall not 
pursue the study of (|2.35|) in this paper. 



From now on, we shall consider the choice Vt a = e Aa ^ s \ which is an hermitian 
and positive semidefmite operator whenever S is hermitian. In this case, H FP takes 
a simple form 



& = S 



— ► H FP =Y, (e Aa(S) - X a ) (2.37) 

aeC 

or in words, H FP is the sum of the free Hamiltonian ( 2.22 ) and an interacting 
Hamiltonian which depends only on the Z-operators. 

The previous results are the starting point of our method. Thus, let us suppose 
we are given a Hamiltonian of the type 

H = -J2X a + U[Z A ] (2.38) 

aeCo 

and whose ground state ip we would like to determine in the usual way: 

Hi* = E ij (2.39) 

A solution of this problem arises if one is able to find an action operator SfZ^] 
which satisfies the Ricatti equation 

U[Z A ]-Eo= J2e A4S) (2-40) 

agCo 



For the eigenvalue problem ( |2.39|) the operator U as defined in eq. (|2.40 ) is a 



datum whereas S is an unknown; however, as it happens in the continuum, one can 
reverse the roles and take eq.( [2.40| ) as defining the interaction term which corre- 
sponds to a given state of the form (|2.24 ). In the sequel, we shall use both kinds of 



interpretations. 

Finally, to complete our generalization of the FP-method we shall study the 
norm of the exponential state (|2.24j) . If S is an hermitian operator, the norm of 
this state can be considered as the partition function of an statistical mechanical 
model with real and positive Bolztmann weights. To show this, we can select an 
orthonormal basis {| /)} of the Hilbert space 7i where the Z-operators are diagonal 

Z A \I)= z AJ | /) (2.41) 
Then the norm of the state ( |2.24| ) becomes: 

Z = (V> | Vo) = (0 | e~ S[ZA | 0) = I (0 I I)\ 2 e- S[ZA ' I] (2.42) 
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In most of the examples (0 | I) will be a constant and that implies that the 
Boltzmann weights are given by e~ s ^ ZA >^. 



Before we apply this construction to particular examples, let us mention that a 
realization of the algebra ( |2.1£| ) in the continuum can be developed in terms of the 
Weyl operators: 

Z? = e m (2.43) 

These operators satisfy the algebra 

X^Z ?2 = e m ^~^Z^X^ (2.44) 

but obviously this expression should be regularized in order to build a consistent 
theory. 

The FP Hamiltonian given by eq. ( |2.3j ) can be obtained just by taking the limit 
h — > of the corresponding fl2.37Q (see Appendix A). Thus our generalization ac- 
counts to a "deformation" of the usual FP-approaches. The word deformation is 
used here with a meaning very similar to the one used in the theory of quantum 



groups JT9[ 



3 The Ising model in a transverse field 



As a prior example of the usefulness of the FP-method outlined in the previous 
section, let us consider the case of the Z)-dimensional Ising model in a transverse 
field (ITF). 

The model is defined on a hypercubic lattice with iV = L D sites and Hamil- 



tonian pl],|22U 



Hlsing,D — Hq + U j 



Ising, D 



a 



where cr?, 



a- are Pauli matrices acting at the site j, with j 
L and \i is any of the lattice vectors \i\ — (1, 0, • • • , 0), ■ 



= Ui 



(3.1) 



]a 



These matrices give a representation of the algebra (|2.18 ) for k 



■,3 d) 

(0,0,- -.,1). 
2. We shall assume 



in 03 



x,z 



for A smaller than a certain D- 



L|) periodic boundary conditions, i.e cx?_ 
It is well known that the ground state of (|3.1[ 
depending critical value A c is non degenerate. This accounts to a disordered regime 
where the order parameter (a?) vanishes. For A > A c , however, the Z 2 -symmetry 
is broken spontaneously and there are two ground states in the termodynamic limit 
with (a-) ^ 0. This is the ordered regime of the ITF model. 

Let us focus our discussion on the disordered region, where the ground state is 
non degenerate. Our aim is, for the moment, to construct this ground state. 

It is easy to check that all the ingredients needed to apply the generalized FP- 
construction are fulfilled by the ITF model. Indeed the algebra of operators X and 



11 



Z are given by the Pauli matrices {crj} and {cr|}, respectively. The vacuum vector 
| 0) is defined by: 



TV 

0) = II I 0), (3.2) 



where 



0?|O>i=|O),., Vj = l,---,iV (3.3) 

The free Hamiltonian ( |2.22| ) is, apart from a constant, the same as the free 
Hamiltonian in eq. ( |3.1|) . The problem is therefore to look for the ground state of 
1]) in terms of the exponential ansatz 

| fa) = e -\SH,-^ N ) | (3 _ 4) 

The corresponding Ricatti equation which fixes S becomes in this context: 

A a i a L + E o = -J2 eAj(5) (3-5) 



where 



A, (S) = l -\S{al ■ ■ ■ , a*, ■ ■ ■ , o* N ) - Sfa, ■ ■ ■ , -o* jt ■ ■ ■ , o* N )\ (3.6) 

The solution of this equation is rather cumbersome, even though in principle it 
can be obtained successively, step by step, as follows. To start with, suppose we 
choose for the functional S in (|3.4j) the expression: 



SF> = -aiY.°toU* ( 3 - 7 ) 

To this choice, it corresponds a state ip^ = e~^ sW | 0) which can be thought of 
as the ground state of a D-dimensional FP-Hamiltonian: 

ttFP,(1) rr , ttFP,(1) _ p FP,(l) 

n Ising,D ~r~ ^ I sing, D ^Isingfi 

N i TV . 

= - 5>J - e Aj(5(1)) ) = i £ #f e^^irf (3.8) 
i=i i=i 

where 

^(i) = ( e -A i( sW)^ _ 1} (3 9) 

and whose interaction potential Uf?£g' D and energy E^^ D are given by: 

N 

vJZ& - e£!& = E e 4 ' 18 "" = E E » 

= (coahaO^Ena-W^C 1 -*! ?^) ( 3 - 10 ) 
j m 



12 



with t\ = tanhai. 

Eq.(piq) yields for the energy E? S ^ D : 

Ef s ™ D = -N (cosh a,f D (3.11) 
and for the Uf s f^ D=12 interaction Hamiltonians: 

N 

uHn%=i = -(coshax) 2 E( 2 ^K+i - ^K+2) (3-12) 

tf£5£U = -(cosh «o 4 E(2ti E - *? E °t°U» 

j A 4 M 

Q+2 2/ Z I Z \ i J-3„Z/ Z Z Z 

+ a i+,i a i-,i a i-^ + Ml ~ M2) + *W+ M1 ^- M1 ^ + A*°J- J ( 3 - 13 ) 



A glance at eqs. (P^) and (|3~T3|) shows that both ^f;^ =1)2 contain a term of 
the form u \ u l+u,- In f & ct, this is a general property of Uf^g D in an y dimension 



D. Then, if we make the identification: 



A = 2(cosh «i) 2D tanh «i (3. 14) 

we may say that Ui S ing,D and uf^g,D coincide to order A but differ in order A 2 and 



higher. On the other hand, if one eliminates a\ in terms of A with the help of (|3.14|) 



and substitutes back the result in ( |3.11| ), one immediately finds the energy of the 
state In the ID case, this energy reads: 

*£2£U = ~ a + v/ttp) (3.i5) 

This expresion agrees to second order in A with the exact result which is given 
in the termodinamic limit iV — > oo by 0: 

e,.,^ = - Jta^i = (1 + A) i; 1; ^) (3.16) 

where F(a, 6; c; x) is the usual hypergeometric function [p4fl . 

From eqs. ( p. 12]) and fl3.13| ) it is clear how to improve the first order aproximation. 



One has to add to the exponent of the ansatz those terms necessary to cancell the 
unwanted ones. In the ID case that corresponds to the choice: 

The parameter a 2 has been introduced as to eliminate the term <T l (T j+2 
appearing m d3~T2l) . Following the same steps as in (|3.10|) we obtain for the potential 



u i%% ) t D=\ and energy E? S ^ D=1 the equation: 



u isin%=i ~ #S%=i = (cosha! cosha 2 ) 2 
13 



x £(1 - ha^ +1 )(l - t x a)a]^){\ - t 2 a*a* +2 )(l - t 2 a*a*_ 2 ) (3.18) 

3=1 



where 



i a = tanha a , a = 1,2 (3.19) 
Now, by demanding Ui S i ng agrees to order A 2 with Uf s f^ D=1 implies: 

E Lin%=i = -iV(coshaiCosha 2 ) 2 

A = 2(cosh«i cosha 2 ) 2 ti(l - t 2 ) (3.20) 
= 2t 2 - t\ 

These equations can be easily solved for Ef s ^g,D=i as a function of A, namely: 

p fp,(2) _ Af A 91 v 

*Isin g ,D=l- 7V tl(A)(2 _ t 2 (A)) V-U) 

where t\(X) is the real and positive solution of the polynomial 

tt + ti + -r- -2 = 0. (3.22) 
A 

In fig. 2 a plot of the values for minus the energy density ef^^ D=1 , ef^g D =i an d 
eising,D=i is presented. Obviously the results improve significantly as we increase 
the order of the aproximation. 

This procedure can be repeated order by order by adding successively terms of 
more complicated structure in the exponential ansatz. The hope is that after a finite 
number of steps the process ends up. This is indeed the case, as can be seen for a 
ID chain with N = 4 (see appendix B). We can say that the ground state solution 
of the Hamiltonian ( |3.1|) can be constructed in a process which is similar to the one 
used in perturbation theory, but with two basic differences; firstly, it incorporates 
at each step non perturbative contributions to the energy and to the ground state; 
and secondly, the process can be completed after a finite number of steps for finite 
chains. 



4 A Fokker-Planck extension of the Ising model 

We have just described an attempt to reconstruct the exact ground state of a 1- 
dimensional Ising model based in the choice of an exponential ansatz. This is cer- 
tainly feasible for small chains, as shown in appendix B, but it becomes an awful 
task for chains already modest in size. On the other hand, it is not our purpose to 
obtain in a complicated manner what can be easily done by other means 0, 0. 
The Ising model in ID is an integrable system and can be solved by performing a 
Jordan- Wigner transformation to fermion operators followed by a Bogoliubov trans- 
formation. In higher dimensions, however, the model is probably not integrable and 
that explains the difficulty in finding an exact result, which quite likely does not 
exist. One of the main interests to seek exact solutions of lattice models consists in 
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the study of their critical properties. Thus, if the model under consideration comes 
from a discretization of a model in the continuum, one would like to recover that 
theory as the critical point of the corresponding lattice model. In any case, the 
choice of a particular hamiltonian is not so important as the universality class to 
which it belongs to. So in this section we propose a FP-hamiltonian which enjoies 
the same universality properties as the 1-dimensional ITF model, while remaining 
close to this latter outside the critical region, in the sense that the energy and other 
physical quantities do not deviate significantly. 

Let us start by considering the most general 1-dimensional ITF hamiltonian 
including all possible couplings involving only two o z matrices, namely: 

N -y N N-l 

Hl sing (Xl, ■ ■ ■ , X N ) = - °j - ~ X X ^rVj<7j +r 

j=l Z j=l r=l 

where, for consistency: 

X r = X N - r , r = l,---,JV-l (4.2) 

Notice that the standard 1-dimensional ITF model corresponds to the choice: 
Xi = Atv_i = A and A r = for r ^ 1 or N — 1. 

We want to find a FP hamiltonian which coincides with ( [4. 1|) except for operators 
involving more than two cr z -matrices, i.e: 



(4.1) 



Ul sing (Xl, • • ■ , X N ) - U Ising (Xi, • • ■ , Aat) C ii,-,i P a 'ii ' ' ' °l v ( 4 - 3 ) 

p>3 ii,—,ip 

According to the FP-philosophy, the knowledge of Hf s f ng is equivalent to the 
knowledge of its ground state, which we shall suppose to be given by: 

| V(«i, ■ • ■ , a N )) = e -* s W'~^) | 0) = e* S-i 1 a ^ +r | o) (4.4) 

where a r = a^-r are assumed to be real and | 0) satisfies Q3.3p. 
The discrete Ricatti equation reads for the ansatz (|4.4| ): 

1 N N-l 

+ 2 E E 

z 3=1 r=l 

N 

-EE cw<--< = -E^ (5) (4.5) 

p>3ii,--,i p j=l 

where 

N-l 

A, (S) = - X a r a*a* +T (4.6) 

r=l 

The cr? operators in eq. (|4.5| ) can in fact be replaced by variables aj taking just 
two values ±1. Then if we multiply both sides of ( |4.5| ) by a convenient number of 
Oj's and subsequently we sum over the 2^ possible choices of a/s we obtain for the 
energy E FP and the coupling constants A r 
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E FP = ~^ E e A ° (5) (4-7) 

{<Tl, — ,<TJV} 



1 



A, = -^ E (Ew^)^ Ao(5) (4-8) 

{<Tl, — .<Tjv} J = l 

Here the subindex in Aq(S) refers to the site "0" (= site "iV") of the chain. It 
is straightforward to prove that 



3 Ao(S) 



N-1 



e~" y "' = J\ cosha r (l — t r a a r ) (4.9) 

r=l 

where t r = tanha r . If we now introduce (14.91) in (14. 71) -(14. 80 and define: 



we obtain: 



1 = -to (4.10) 



_J^FP N-1 

z F N p = ——= J] cosha r (4.11) 



N-1 



A r 

e N n=0 



EWr-n , (r = !,■■■, N-1) (4.12) 



where e N p is minus the FP-energy density of a chain of length N. 

The structure of eqs.( [4.1lD -( |4.12|) has already emerged in, for example, eq. (|3.20D 



(see also eqs. ( |B.13| )-( [B.14| ) of appendix B). These equations relate the t n (or the 
a n ) to the coupling constants A r and, therefore, one could in principle relate the 
FP-energy E FP with the A's. In this whole process, the parameters Ci u ... t i which 
stablish the difference between our model and the pure Ising model ( |4.1|) play no 
role. In fact, they can be computed a posteriori in terms of the a's. This means 
that our prescription to associate a FP-model to the Ising is well defined. 

The next step is the resolution of eqs. fl4.11|) - fl4.12|) . For that, let us introduce 
the Fourier transformed t n of the t n : 

i N-1 

*» = -nr T E e^tm (4.13) 
ViV m=0 



and consequently eq.( |4.12j ) becomes: 



l^ = -Ee~(U 2 (4-14) 

e N m=0 



Then, if we apply an inverse Fourier transformation to (14.140 we get the £ n 's in 
terms of the A's and a parameter Cn, namely: 

i N-i \ 

L = -^=e m (C N - E ^p) 1/2 (4-15) 

where e m = e^-m = ±1 and 
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N-l N-l 
n=0 



m=0 



is a real constant which can be fixed by imposing the eq. (|4.10| ): 

N-l 1 N-l N-l 



1 1 

= ^ im = TT ^ ^m(CiV ~ e ~ 

N±?n N ^ 



A, 



vl/2 



(4.16) 



(4.17) 



m=0 J ' m=0 r=l TV 

One can easily realize that the trivial solution for the case A r = 0, Vr, corre- 
sponds to: 



C N = l 



1 



N 



e m = l (4.18) 

The presence of interaction modifies these values; for convenience, however, we 
shall continue maintaining e m = 1. 

When A = Ai = Xn-i, (K — 0, otherwise), one recovers the usual ITF model 
4]). Eqs. (|4.13| ) and ( |4.17|) becomes now: 



1 N-l 

J. x — > 2 7T i n m 



e - (ON--?? 

iv m=0 e N 



2A .2irm 1/2 
cosf— —)) 1 



N 



1 



N-l 



iv m=0 



2A ,2nm 1/2 
cosf^^)) 1 ^ 



iV 



(4.19) 
(4.20) 



Notice that ( 4.20Q is nothing but the n = case of (|4.19 ). 

We shall look for solutions of ( |4.19| )and Q4.20| ) for which the t n are real quantities. 
This is guaranted if we impose the constraint: 



< u N 



2A 



FP 



Cn&n 



< 1 



(4.21) 



The solution of ( |4.19| )and ( |4.20| ) involves, in principle, four steps: i) to express 
Cjv as a function of A and e^- p , ii) to rewrite t n in terms of A and e^- p , iii) to introduce 
the t n in eq. Q4.lip as to obtain an equation for e N p as a function of A, and iv) to 
solve this equation. 

Fortunately, this rather complicated elimination process can be explicitely ac- 
complished in two extreme cases according as: a) un = 1 and N arbitrary or b) 
un < 1 but iV — > oo. Next we shall describe our results leaving the details and 
derivations to the Appendix C. 



a) Case — 1 

The values of t n and eff are given by: 



(sin 



7T ^2 

27V ) 



sinf (n + i)sin^(n - \) 



D FP 

'TV 



V2 



N sin sin 



(4.22) 
(4.23) 
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where the angle On in (|4.23|) is defined as: 



cos^jv = 2cos(^) - 1 (4.24) 

If we take the limit iV — > oo in eqs. (|4.22|) and ( [4.23| ) we deduce that the model 
exhibits critical behaviour. First of all, from eq . ([4 . 2 2|) we obtain 

lim t n = (4.25) 

JV^oo An 2 - 1 

and hence for n large the "amplitudes" t n have a scaling behaviour. It is curious to 
observe that the value of t n given by eq. (|4.25|) coincides with the correlator of two 
a x, s matrices in the standard ITF model at the critical point H: 



We may wonder whether these two results are related. 

Another indication of criticality arises from the 1/N 2 correction to the ground 
state energy density, which is given by: 

^-~ P |l + ^(l-^co t (^)) + 0(] l)] (4.27) 

where 

p P 2V2 3.554 

C P = . / . m = (4-28) 

7rsm(7r/v2) 71 

Recall that e N p is minus the energy density of the ground state. 
It is well known that in a ID conformal invariant theory this correction satisfies 
the Cardy formula: 

KCV S 

e = e OQ + — (4.29) 

where v s is defined through the dispersion relation of the elementary excitation 
above the ground state, i.e. u(k) ~ v s k, and c is the central extension of the 
Virasoro algebra underlying the corresponding conformal field theory. 
Compairing eqs. (|4.27|) and (|4.29|) we learn 

c^ = ^(l--^cot(^)) (4.30) 

Hence we would need to know the value of v s in order to derive the central 
extension c (see section 5 for some considerations concerning this point). Our ex- 
pectations are that c = ~, in which case we would be actually describing an Ising 
model; but we shall arrive to this result by studying other quantities, such as the 
critical exponents. 

It is amusing to observe that the value of the critical energy density ( |4.28| ) is 
comparable with the one of the ITF model: 
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JTF 



4 

7T 



(4.31] 



To complete this comparison, let us recall that in ID the critical value of A is 



given by A 



ITF 

c 



1, whereas in our FP-ITF model we find 



A 



FP 



71 



4 v / 2sin(7r/v / 2) 



0.6979 



(4.32) 



b) Case Uoo < 1 

From the discussion of the previous paragraph we expect that this case will 
account for the non critical regime of the model. The expressions for the t n and the 
energy e^f are given here by (see Appendix C for computation details ): 



7T 



2 5/2 



2A, 



11 + Ah) 2 



r 



■n 1 



im + ±) F(n-\, n+i 2n+l; 



4A F 



r(n + i) 



;i + a b ) f(-1,1,i 



ill. 4A fl 
2' 2' X > (1+A H ) 2 



4A 



l 2 



2' 2' '(1 + A^ 



where we have introduced a new variable A^ through: 



2X, 



Ur 



1 + A 2 



(4.33) 
(4.34) 

(4.35) 



R 



Equations ( fl.33| ) and ( 4.34[ ) make clear that \r rather than A is the effective 
coupling constant which governs the behaviour of the ground state. In a certain 
sense, A# results from the dressing of the "bare" coupling constant A by all the 
interactions present in the FP-Hamiltonian. The explicit form of A as a function of 
A# is very complicated and it will not be needed; in any case it can be obtained by 
comparing eq.( [4.35| ) with the eq.( |C18| ) of Appendix C for e^f as a function of the 
t„'s. In fact, it is not difficult to show that 



A = A fi (l--A| + 0(A 4 /? )) 



(4.36) 



and 




(4.37) 



Making use of (|4.37|) one sees that is not singular at the critical value 



X R — 1. However, the "specific heat" defined as the second derivative 
actually exhibit a logarithm divergence 

d 2 e^ SA, 



d 2 t 



does 



V„FP 

~d\7 



7T 



log(l - \r) 



(4.38) 
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This is precisely the same logarithm singular behaviour that one would expect 
from the original ITF model ||. This implies that the critical exponent a defined 
through ~ (1 — \) a is the same for the two models, i.e. a = 0. 

Let us concentrate now on the behaviour of the t n for n large. In the critical 
case Xr = 1, we have seen that they satisfy the scaling law (|4.25|) ; if A# < 1 we 
obtain instead 

£ - ^72 A S for n » 1, X R <1 (4.39) 



and so eq. (|4.39| ) allows one to define a correlation length £(A#) through t n ~ e <= 



as 



r 1 (A /? ) = -logA fi ~(l-A i? ), A R ~1 (4.40) 

which implies in particular that the critical exponent v is equal to 1, just like in 
the standard ITF model. This is again another indication that the model we are 
considering belongs to the same universality class as the standard ITF model. 



5 A change of variables: excited states 



Until now we have considered the construction of the ground state of a lattice Hamil- 
tonian in the spirit of the FP procedure; at this point, it seems natural to wonder 
what can tell us the method about the rest of the eigenstates of the Hamiltonian. 
This is, of course, important both for the study of gaps in the spectrum as for the 
dispersion relation of elementary excitations. 

In order to fix ideas, let us start by posing the problem in the context of Quantum 
Mechanics. Hitherto, we have been using the fact that a wave function of the form 



-S(x) 



x G R, is a ground state of the FP-Hamiltonian: 



H 



FP 



d 2 



dx ' 



cPS_ 

dx 2 



(5.1) 



If, for example, we take S(x) 



1 r 2 

2 ' 

{R FP 



then H FP coincides with the number 
operator for the harmonic oscillator (H i ' F = a^a = Af, [a, a'] = 1). If we want 
to study the excited states of Q5.1J ), we could make the ansatz: ip(x) = A(x)e~ s<yX \ 
which can be thought of as a change of variables from ifj(x) to A(x) since the function 
S(x) is known a priori. The eigenvalue problem for if)(x) becomes then: 



H FP i)(x) = AEip(x) 



where 



H RFP A(x) = AE A(x) 



(5.2) 



H 



RFP 



d 2 | 2 ds d 

doc doc doc 



(5.3) 



and AE represents the difference of energies between the excited state ip(x) and the 
ground state ipo(x) -whose energy is actually zero. We shall call Hrfp the reduced 
Fokker-Planck Hamiltonian. The ground state solution using the new variable A(x) 
corresponds to A(x) = Cte. In particular, for the harmonic oscillator, eq.( |5.2|) reads: 
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1 d 2 A dA . _ . , w lS 

~ 2*? + < 5 ' 4) 

which becomes a confluent hypergeometric equation after the substitution y = x 2 , 
and whose normalized solutions are given by the well-known Hermite polynomials. 
It is worth noticing that H RFP is not an hermitean operator with respect to the 
measure /i = dx, but it is so for the measure fi = e~ 2S ^dx. 

Next, we shall try to generalize the previous methods to more complicated Hamil- 
tonians. More precisely, the problem we face is the construction of the spectrum of 
the FP-Hamiltonian associated to fl4.1|), namely: 



N N N-l 

H FP = - ZK - eAj(5) ] = " £[°i - ^ II (1 - *r^+r)l (5-5) 
j=l j=l r =l 

whose ground state ipo was defined in 

I N-l N 

I V^o) = exp(- J2 E (*r(TjCr* +r ) \ 0) (5.6) 

4 r=l j = l 

In ( ^.5|) e^- p represents minus the FP-energy density given by (|4.11|) and t r = 
tanha r . In principle, this ground state may not be unique. For example, this 
would happen if the Hamiltonian ( |5.5|) enters into an ordered region characterized 
by ground states with opposite non-vanishing magnetization. We shall not discuss 
this possibility here. 

The change of variables adequate for our purposes corresponds to: 

\^)=A(a z ir --,a z N )\iP ) (5.7) 

where A(af, • • • , <j z n ) is the most general operator involving a z, s. With the notations 
of appendix B: 

A(a z 1 ...,a z N ) = J2A Ia a z Ia (5.8) 

a=l 

where in contrast to the appendix, we shall allow in I a the presence of any number 
of l's, independently of their parity. The eigenvalue equation (|5.2|) reads now: 



H RFP | A) = AE | A) (5.9) 

where 



A) = A{al,...,a%)\Q) = Y,A Ia \I a ) (5.10) 



a=l 



and 



N N N-l 

H Rfp = Ee A, (s)(1 _ a?) = II (1 " " (5-11) 

j=\ j=l r=l 
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For reasons will be evident latter, it will be of some convenience for us to intro- 
duce the parameter: 

AE 

e =2^ (5-12) 

and in what follows we will also assume that H RFP is divided by 2ej^ p to give H RFP . 
Notice that in the absence of interaction (i.e t r = 0, Vr), we have: 

fyWRFP = J2^—^=M (5.13) 

3=1 2 

and so H° RFP coincides with the number operator A/", i.e, the operator which in 
this context counts the number ni a of l's of a given state | I a ): 

M\I a )=n Ia \h) (5.14) 
An interesting property of the operator H RFP is 

2^ N-l 

Y,(Ia \ H RFP \ h) = U(l-tr)n Ih (5.15) 

a=l r=l 

This equation suggests some kind of connection between the evolution described 
by the Hamiltonian H RFP and a Markov process. As a consequence of (|5.15 ) we 
find 



JV_i TIT At 

e=U(l-t r ) ^ f ( 5 - 16 ) 

r=l l^a=l A Ia 



A useful way of reexpressing H RFP is 

N-l 

({^r}) ^ t t r A. r -\- ^ tr\tr2^-Ti,T2 ' ' ' 

r=0 l<ri<r-2<iV— 1 

• • • + ( 1) ^ ] t ri • • • t ri A rij ... jr; -(-••• 

l<r 1 <r 2 <---ri<N-l 

- tit 2 ■ ■ ■ t N -t A 1)2 ,...,n-i (5.17) 
where the operators A's are defined as 

N 1 - a x 

3=1 Z 



Notice that the first sum in ( |5.17| ) starts at r = and so reproduce the operator 
M (recall that t = -1). 

The action of A on a generic state | I a ) —\ • • • , fi N ), (/^ = or l,mod(2)) is 
given by (see fig. 3): 



N 

A rx ,-, n | Hi, ■ ■ ■ , Hn) = bi>H I ' ' ' ' to + ^ ' ' ' ' + 1 ' ' ' ' ' + ' ' ') ( 5 - 19 ) 

7 = 1 
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The diagonalization of the Hamiltonian ( |5.17[ ) seems, a priori, not much easier 
than the diagonalization of the Hamiltonian (|5.5|) . However, we have developed a 
real renormalization group approach which enables us to deal with this apparent 
difficulty; but the results concerning this point will be presented elsewhere PB] . For 



the moment, and with the purpose of getting some insight into the Hamiltonian 
( |5.17| ), let us consider only the first term in it which is linear in the i's, namely: 



N-l 

ffW Rfp = _j2 tr A r (5.20) 

r=0 



It is easy to see then that the state: 

N 



3 * 

where | j) =| ■ • • 1 • ■ • 0) is an eigenstate of RFP ; 



H w RFP | V™> = -VN t m I Vm> (5.22) 
Hence, eq. (|4.15| ) can be used to obtain the dispersion relation for a plane wave 



of momentum k = 

AE FP (k) = 2e F N p (c - V2 (5.23) 

The gap is obtained for k = and vanishes in the critical point, i.e when = 1. 
Moreover, for k small the dispersion relation (|5.23| ) is linear in k: 



AE FP (k) ~ \ e F Jk (5.24) 

whence we deduce a speed of sound = | e^ p . This confirms our expectations 
that the value un = 1 actually corresponds to a critical point. Unfortunately, the 
speed of sound we derive from eq.( |5.24j) is not the speed of sound of the excitations 
of the whole Hamiltonian (|5.17|) . In fact if we take for granted that c = 1/2 is 
the value of the central extension of this model, then the speed of sound for the 
Hamiltonian (|5.20 ) should be given by: 



^^(1-^cotQ) (5.25) 

This equation suggests that the role of the additional terms in the Hamiltonian 
( |5.17|) is simply to redefine the value of Vf\ 



In section 4 we introduced an effective coupling constant A# in terms of which the 
energy e pp displayed a critical behaviour. If we now rewrite eq. (|5.23|) as a function 
of A r we obtain 



AE FP (k) = 2 J (1 + A| - 2X R cos k) l/2 (5.26) 
V Aij 

It is noteworthy to compare this relation with the exact dispersion relation of 
the standard ITF model. As it is shown in [Rfl 
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AE e j^(k) = 2(1 + A 2 - 2 A cos k) 1 ^ 2 (5.27) 

The evident similarity between eqs.( |5.26| ) and (|5.27| ) is again another confirma- 
tion that the model we are describing is very close to the usual Ising model; indeed, 
we believe they are identical at the critical point. 



6 Conclusions and final remarks 

This paper represents a first step in the application of the Fokker-Planck techniques 
to the study of quantum lattice Hamiltonians. The main idea is to approximate 
the Hamiltonian of interest if by a FP-Hamiltonian H FP which shares with H a 
certain number of properties. Firstly, we have considered FP-Hamiltonians which 
agree with H to first and second order in perturbation theory. Next, and more 
interesting, we have built up a FP-Hamiltonian belonging to the same universality 
class as H. In both cases the advantage of working with a FP-Hamiltonian is that 
we know the exact ground state of this latter. In the example we have considered, 
i.e. the ITF model, we have been able to solve the equations fixing the ground state 
parameters in terms of the coupling constant, and that makes feasible the study of 
the critical behaviour of the model. 

Below we list some issues we believe deserve further investigation. 



Modelistic: 

• The Ising model for D > 1. Sections 4, 5 and 6 can automatically be extended 
to ITF models in higher dimensions. As a matter of fact, increasing the space 
dimension does not represent a priory a barrier to the FP-techniques, as is 
shown in appendix D. Our feeling is, however, that realistic results would 
require of exponential ansatz involving more than two o z operators. In this 
sense, the analysis of the D = 2 ITF model still remains as a challenge in view 
of its connection with the 3D classical Ising model. 

• Zn- gauge models in D = 2 and 3. Here we can inmediately apply the FP- 
procedure developed in the paper to construct the ground state in the confining 
phase; in 2D this phase is dual to the ordered phase of the 2D ITF model p6|. 



The Potts model. FP general machinery of sections 2 can be easily applied to 
the case of the ID Potts model . The Potts model is a fruitful generaliza- 



tion of the Ising model and has been received much attention in the context 
of statistical mechanics. In the notation of section 2, the k-state ID Potts 
Hamiltonian is defined by: 

HpottsW — Ho + \Up ot ts(Zi, • • • , Zn) 

k-l N k-l N 

= -EE^;-^E(¥j + i) n (6.i) 

n=l j=l n=l j=l 
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where the operators X and Z verify XZ = e^ZX and XZ^ = e~^~ Z^X. 
The Hamiltonian (|6.1|) constitutes an example of a traslational invariant 
and hermitic operator, which includes the Ising Hamiltonian as a particular 
case (k = 2). 

Fermion models and Heisenberg models. In the example considered in this 
paper, the operators Vj entering in the exponential of the ansatz commuted 
among themselves. This property facilitated a lot the computations. However, 
this is not the case for models with fermions or spins P7J. Hence, the question 
of a convenient implementation of the FP-techiques so as to include fermions 
still remains to be done. 

Degenerate ground states. The general FP-construction introduced in section 
2 applies also when the unperturbed hamiltonian Hq has degenerated ground 
states, as is the case of the ITF model in the ordered phase; here the compu- 
tation goes in parallel to the well-known degenerated perturbation theory of 
Quantum Mechanics. 

Universality: 

Apart from the critical exponent a and v computed in this paper, one would 
like to determine the rest of exponents, specially the fractional ones. In fact, 
only after this determination had been accomplished, could a complete agree- 
ment with the critical Ising model be claimed. 

If our statement concerning the universality class of the FP-Hamiltonians is 
correct, then we have a method to compute the exact values of critical expo- 
nents. One may wonder how the scenario would be like in dimensions higher 
than one, where non trivial integrable models do not exist. 

FP-procedure seems to throw a light on the way we could compute the corre- 
lation functions of some operators. For that purpose, we could use the ground 
state energy as a generating function of correlations along the lines of the 
Hellmann-Feynman theorem. 

A fundamental question concerning the FP-aproach undertaken in this paper 
is how to find a "minimal FP-ansatz" capable of reproducing the same critical 
behaviour of the original Hamiltonian H. We do not have a neat answer to 
this question yet; however common sense recommends to begin with "simple" 
ansatzs and from these to construct more complicated ones. The final objec- 
tive must be, of course, to capture the essential physics of the model under 
consideration. 

Integrability: 

An interesting open question is whether the FP-Hamiltonians are integrable 
or not, since this would explain some of the magic these Hamiltonians enjoy. 
If so, one would like to know whether or not there exists a kind of Yang-Baxter 
equation underlying or responsable for such a behaviour. 
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A A q- deformation of the FP-quantum mechan- 
ical Hamiltonians 



In Section 2 we were concerned with a Fokker-Planck approach to Hamiltonians 
which appear in continuous or lattice Quantum Field Theories. We devote this 
appendix to show that most of the formulae appearing there have an analogue in 
ordinary Quantum Mechanics. For pedagogical reasons, let us concentrate on ID. 

To begin with, let us write down the analogous to eqs. (|2.3|) , (|2.4| ) and (|2.5|) for 
a wave function ip(x) = e~ s ^ x \ xGR: 

jjFp _ 1 r d 2 ^ ,dS. 2 _ d 2 S .. 
2 dx 2 dx dx 2 

= \R ] R (A.l) 

R = i-^+if- (A.2) 
dx dx 

Then if we put H FP = H + U FP (x) — E FP , we obtain the Ricatti equation 



U FP (x)-E FP = -U^) 2 -^] (A.3) 



1 . dS .2 d S ~, 

Now suppose that x = 9 is compactified on a circle of radius 1. Then it makes 
sense to work with the variables Z = e ld and rewrite the wave function ip(9) as a 
single valued function of Z, i.e ip(Z). Similarly, instead of working with the usual 
momentum p = —ih-^, we may choose the operator X = e ip . Both X and Z satisfy 
commutation relations 

XZ = qZX , q = e ih (A.4) 
Now we define the so-called g-derivative D q 

„ X — X ^ , . . 

D q = i A.5 

q-q 1 

which coincides with the usual 4- in the limit h — > 0. The free Hamiltonian 

dx 

H = can therefore be replaced by a "g-deformed version": 

Ho(l) = -\Dl = 2{q ^ q - 1)2 (X 2 + - 2) (A.6) 
and the corresponding FP-Hamiltonian becomes: 

H FP ( q ) = _ X [X 2 + X- 2 - e^)- 5 ^) - e^-s^z^ 



\R\q)R{q) (A.7) 



with 
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R(q) 



q-q 



_^[S(q 2 Z)-S(Z)] X 2 _ e -l[S(q 2 Z)-S(Z)}^ 



(A. 



where the hermiticity of S has implicitely been assumed. 

Finally, if we write H FP (q) as H (q) + U FP (Z) - E FP , we get 



U FP (Z) - E 



FP 



1 



■[2 



D S{Z)-S{q 2 Z) _ p S(Z)-S(q- 2 Z)l 



(A.9) 



2{q-q-^r 

Notice that this Ricatti equation determines the value of S at the point q 2 Z as 
long as the values of 5* at the points Z and q~ 2 Z are known. In this sense, ( |A.9| ) 
should be considered as a second order finite difference equation. Besides, either 
( |A~7D , flOp and ( |A~9| ) become (|ATT|) , (TA~^ ) and (TAT3| ), respectively, in the "classical" 
limit g — > 1. Moreover, if g were a /c-th root of unity, then the 9 variable could be 
restricted to take a discrete set of values, let us say 6 = , n = 0, 1, • • • , k — 1, 
in which case we would obtain a Z^-invariant theory. 



B Exact FP-computation of D=l finite chains 

A suggestive feature of the FP-approach considered in this paper is that it provides 
a recipe to compute exactly the ground state and energy of a whole family of 1- 
dimensional hamiltonians defined on a finite lattice, among which the ITF model of 
section 2 would constitute a particular case. Here we present the derivation of this 
result. 

Firstly, let us consider the lattice configuration as a collection of 0's and l's. 
The presence of a 1 (respec. a 0) at the site j of the lattice should be thought of as 
equivalent to the existence (respec. absence) of a a z at this position. Besides, let 
us assume that the number N of sites of the lattice and the number of l's arranged 
on it are both even. Under these assumptions, it is not difficult to convince oneself 
that there exist just 2 JV_1 different possibilities to distribute a set of 0's and l's on 
the lattice. We can associate to each one of these possibilities a label I a , namely: 



{I a f a=l = {h = (0, 0,^ • } ; 0), I 2 = (1, 1, i • 0), • • • , I^-x = (1, 1, ■'; 1, 1)} (B.l) 

Indeed, we can think in the labels I a as the elements of a certain equivalence 

class -|=— . Here Z2 reflects the sole presence of 0's and l's in the lattice, with the 
usual rule of sum: 

1 + = 1 mod 2 , 1 + 1=0 mod 2 

whereas jF a takes into account the fact that the number of 1 is constrained to be 
even: 

N 

J fl = ^/ a (j) = 0mod2 

i=i 

Secondly, to each label I a we attach an operator a z Ia defined as follows: 
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N 



a 



l = U(Pi) Ia<i) ( R 2) 

3=1 

( to I a — (1,1,0,0,1,1) it would correspond the operator a z Ia = afa^cr^aQ, for 
example). 

It is easy to see that the local operators o\ a verify the fusion rule: 

(B-3) 

and, of course: 

K) 2 = 1 (B.4) 
Likewise, the following commutation relation is satisfied: 

a*al = (-1)^(3)^ (B.5) 

Finally, let us introduce a kind of generalized ITF hamiltonian, characterized by 
the fact that the interacting term includes couplings to all scales and to any number 
of spins: 

H = -J2^j - E Aj B of B (B.6) 

j=l a=2 



The question of finding the ground state and energy of the hamiltonian (|B.6| ) 
can be faced as follows. Suppose we assume for the first the ansatz: 



but now with: 



a=2 

and where the state | 0) verify (^.3|). 

We want to obtain the FP hamiltonian associated to the ansatz (P-7| ), so we 
apply H = — Z)fLi of on it; as to get: 

N N ^N-l 

R FP = _ £ ((J . _ e A,( S)) = _ _ e -E a . 2 UJM a) (B g) 

3=1 3=1 

A lengthy but straightforward calculation, where the use of eqs.( B.2| )- (|B.5|) is 
essential, allows one to rewrite ( |B.9| ) as: 



^ P = -Ekl-^E<{E^ E (b.io) 

a=2 ^ P ' a u ...,a p = 2 

-'ai "I - ' ' ' I<ip -*a 

where we have introduced 
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A 3 = II COsh [ J a(j>/J 



a=2 



t { p a = Ia{j)ti a = I a (j) tanha 7a 



(B.ll) 



(B.12) 



In fact, if we take translational invariance for granted, the value of Aj remains 
invariable whatever be the site j on the lattice; we may then express A,- = A, Vj = 
1, • • • , N. Assuming that this condition is fulfilled, the energy associated with ( [B.10|) 
is given by: 



A' 



E 



FP 



j=i P >i y- 



1)p 



E 



tf.-.tf- 



(B.13) 



I<n H " + Ia p — h 



Moreover, the comparison between the H FP of (|B.10|) and the generic hamilto- 
nian ( |B.6| ) suggests the identification: 



N 

i=i P >i 



-i)p 



E 



tf-~tf , a = 2,---,2 N - 1 



(B.14) 



To make explicit these formulae, let us apply them to the ITF model of section 
2, but now for the case of a finite lattice with N = 4 sites. In table 1, we have 
summarized the notations and coupling constants appearing in this CcLS6, clS well as 
the interaction terms these constants give account of. 

A glance to table 1 suggests the choice 

| = e l"l E'=l CT M + l + 3°2 Ei=l <^J +2 + ^3^<7§<Tfa! | ^ (B.15) 

where the factor | in (P-15| ) has been introduced in order to avoid overcounting in 
the corresponding term (for N = 4, Y,j=i ffj&j+i = Ej=i o"|cr|+3)- 
Now we require that Hf p F coincides with Hitf- 

Hf T p F = H ITF = - E a i - A E ( B -!6) 

3=1 3=1 

but since we have three different couplings it must be 

ti 2 — ti z = = th = ti = tanh«! 

tie = ti 7 =t2 = tanha 2 (B.17) 
ti 8 = t 3 = tanh «3 
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On the other hand, one is free to select a particular site, say j = 1, to 

1 



ret: 



A = A l = \{ cosh[/ a (l)a/J 



(l-*I)(l-*I)ni-*§) 1/2 
Now, introducing eqs.( |BT7D - (pig) in (gig )- (gig) we obtain: 



rpFP 
ITF 



1 + tlt 2 t 3 



and 



A/ 2 — Aj 3 — Aj 4 — A/ 6 



4 (l-tf)(l-tl)V2(l_ t |)l/2 

2ti(l-t 2 )(l-t 3 ) 



A 



\j a — A 



(l-t?)(l-t|)V2(l_ ti) l/ 2 

o = (i-t 3 )(t?-t 2 ) 



A/ 8 = = t 3 + £ 2 t 2 



Eqs.( |B.19] )- (|B.20D can be easily solved, and in fact we obtain: 

FP _ (1+^I) 1/2 

°ITF — 1 7 
1 — t 2 

where t 2 is the real and positive solution of the polinomial: 



'-0 



^-2^-^| + l = 



namely: 



1 



:[A 2 + 1 + VTT¥ - 2(A 2 + v/TTA 1 ) 172 ] (B.23) 



t 2 (\) = 

Now if we substitute eq.( |B.23j ) in ( |B.21[ ) we find the final result: 



(B.18) 



(B.19) 



(B.20) 



(B.21) 



(B.22) 



effpW = i[ v / l + A 2 + Av / 2 + \/l + A 2 -AV2] (B.24) 
This expression coincides precisely with the exact value given by the formula : 

(B.25) 



ciTF(A) = i fc ^l + A» + 2Aoos[S^±^] 



when we particularized to the lattice with N = 4. 

C Explicit computations in the FP- version of the 
ITF model 

In this appendix we deduce the formulae for the t n and ef/ 3 exhibited in section 4. 

The case ujy = 1 is very easy to handle because the square roots appearing in 
eqs. Q4.19D - (|4.20D can be performed in terms of trigonometric functions. In particular, 
we get from eq. (|4.15|) 
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t. 



2Cn . nm, 
~N~ Sm( AT ' 



and introducing this equation into eq.( j4.17|) we find: 

N-l 



1 



V2C N W . .Tim. 
> sm 



N 



m=0 



N K 2N' 



(C.l) 



(C.2) 



From this euation the value for Cn follows 



C N = \{N\m^Y (C.3) 
On the other hand, working out eq.( |4.19D with the help of ( |C.3|) and ( |C.1| ) yields 



t n (N) 



sm 



2 ( TT - 

• 2N . 



sin[f + sin[f(n + |)] 
In order to compute the energy efj P we rewrite eq. (|4.11|) in the form 

N-l 

4 P = IK 1 + 

n=l 

Now we can make use of eq. (|C.4|) to express the monomials of (|C.5|) as 
1 + t n 



sm 2 (f ; 



1-U 



sin [f (n + I)] sin [f (n + 1)] 



sin (f + ^)sin (f - e -f) 



n sin[f (n+ |)] sin[f (n + |)] 

where 0jv was defined in equation ( |4.24| ). 
Now using the formula 



N-l 



-p-p 7T 1 sin(7ra 
Y[ sm ( n _ a )= 

n=l JV z 



sin(f: 



we obtain 



iV-l 



II(l + tn) = (iVsin -?-) 2 

n=l /iV 



JV-1 



IK 1 -*" 



/sin ^ sin^ 



n=l 



sm 



From these equations and eq. Q4.23|) the eq.(|(J.3|) follows directly. 



(C.4) 



(C.5) 

(C.6) 
(C.7) 



(C.8) 

(C.9) 
(CIO) 



Next, let us consider the case where N — > oo and < < 1. Now the sums 
can be converted into integrals, as for example eq. Q4.19 ): 

*■=-/. i )1/2 
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^ n+l VcZ r dx cos ( nx )( 1 + MooCOSX )i/2 ( C11 ) 

7T JO 



But replacing cosx = 1 — 2sin (|) in ( |C11| ) and subsequently changing vari- 

X 

2 



ables | — > x we get: 



-i 2 /TwTT - T f* „„„/<-> — \ /1 2u °° ,.;„•■> ,.j ^ 



t n = (-l^-VCooa +Woo) / cos(2nx) (1 - --^sin^) 1 / 2 (C.12) 

7T v JO 1 + Mqo 

This integral can be performed if we first expand the square root according to 
the formula: 

(i -zr = t w r( T r rtn zm ( c - 13 ) 

r(-a)r(m + 1) 

and then we use the result 

„2L r (-1)" r(2m+l) ^ 

/ 2 rfx cos(2na;) sin 2m x = I 22m r(m+n+i) r(m-n+i) m - n ( C14 ) 
Jo [0 m < n 

Finally, we obtain the expression 



The equation which yields is obtained by setting to = — 1 in eq.( |C.15| ), 
namely 

CZ(1 + Moo ) 1/2 F(~ i, 1; = 1 (C.16) 



Taking the square of eq.( |C.16|) and using the fact that = — , we deduce 

t?oo Woo 

an equation for e^fas a function of A and u^: 



FP n\(^~ JrU °°\ T?2l ^ - 



«£f = 2A(— ^) F'(— , -, 1; — (C.17) 

Replacing now in terms of the parameter A# introduced in ( |4.35| ) one obtains 
eq « . 

As for the equation ( |4.33[ ) for the £ n 's, it is easily derived from ( |4.35| ) and ( |U.15p - 



( U.16|) . Likewise, expression ( l.34|) is obtained by taking the N — ► 00 limit of 



my- 

00 

i^na-o- 1 (cis) 



n=l 



and using (|C.15| ). Notice that the exponent in eq.( |C.18| ) is —1 and not — ~ as in 
eq.( |C5| ). This takes care of the constraint t n = t^-n which in the N — > 00 limit 
becomes t n = t_ n , with nGZ. 
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Eq. flCTTSP for 



1 gives the value t r 



If Mno < 1 we can also find 



4n 2 -l' 

an indication of the behaviour of the t n , just on condition that n be large enough. 
To this end, we can carry out a saddle point evaluation for t n on the basis of the 
standard integral representation of hypergeometric functions [p4" |. Indeed, we can 
write ( |C.15|) as: 



Coofl +U C 



2tt 



\n-2) 



where we have introduced the parameter 



f 1 dt 


[*(!-*)] 


10 


_l-zt 



(C.19) 



2w r 



1 + u c 



(C.20) 



Due to the particularly compact form of the integrand, the saddle point in 
eq.( U.19 ) turns out to be independent of n. Our result is 



t ~ 



1 c n 



1 + Jl 



x 



ut 



1 + M00 + Wi 



1 + M c 



1 + ^1 



n»l (C.21) 



Finally, the expression ( [4.35|) can be used to replace Uoo by \r and in this way 
we get eq.( [4.39| ). 



D Higher dimensional ansatzs 

One of the reasons to justify the introduction of new ansatzs for the ground state 
of the ID ITF model, amounts to our wish of discussing higher dimensional models 
for which exact solutions are not known. In this sense, the obvious generalization 
of the l.D-dimensional ansatz proposed in eq.( f4.4|) becomes 



|^R})> = ex P (i££a r a?a? +r ) | 0) (D.l) 

V reA'jSA / 

where A is the hypercubic lattice and A' = A — {0}. We shall suppose that 
a r = a_ r with r G A'. 

The FP-Hamiltonian corresponding to ( p.l[ ) has the form 

H FP = _ Y^j* _ g-ErgA'* (D.2) 

jeA 

Expanding eq.( p.2| ) we arrive to a generalized ITF model 

H FP = _ E FP _ £ a , _ 1 ^ ^ + + . . . (D _ 3) 



jeA 2 j G A 



with 
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E 



FP 



-Ne 



FP 

'N 



FP 

N 



-N Y[ cosha r 



reA' 



(D.4) 
(D.5) 



neA 



and 



to = -1 



t r = tanh a r , r ^ 



To solve eqs.( D.4 )-( D.5 ) we introduce the Fourier transformation 

1 



tr 



E 

meA* 



„27ri ■ 



~ t 



(D.6) 



(D.7) 



where A* is the dual lattice of A. For a hypercubic lattice both A* and A coincide. 
Following the same steps of section 4, we find 



L D / 2 



reA' 



A„ 



1/2 



-A 



where 



C 



N 



(D.8) 



(D.9) 



meA* neA 

and e m = e_ m = ±1. For convenience, we shall fix e m — 1- 
In the limit N — ► the t n are given by 



2^ djfej 

o 27 



2^ dA; 



o 2tt 



O gink 



Coo 



>FP 



K cos /c a 



^oo a=l 



1/2 



;d.io) 



where A a is the value of the coupling constant multiplying the operator cr|c"| +Ata in 
eq.( |L).3| ) (fi a is the unit vector along the a t/l -direction). 

In the special case of D = 2 the integrals of (p.lO|) can be explicitly determined 
in a way which goes parallel to the followed in appendix C. Writting n = (711,712) 
and making use of eq.( p.6|) we obtain after a lengthy computation 

2(Ax + A 2 ) 



x 



^1 + Aij x + \r 2 ) F 2 f— — , -, -, 1. 1 



^(A^ + Xr 2 ) 

4A fil 



4A 



P-2 



2' 2' 2' (1 + A Bl + A K2 )2' (1 + A Rl + Ar 2 ) 2 ' 



D.ll} 



and 
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4A 



,1 + Aitj + A_r 2 ^ 



t 



1 

""" 2 ~ 2V? 

2 r(n x + n 2 - 1) r(m + 1) r(n 2 + 1) 



v l + A^ + A_r 2 " 



r(2m + 1) r(2n 2 + 1) 



F 2 (ni + n 2 — |, ni + §, n 2 + |, 2rii + 1, 2n 2 + 1, 



4A F 



4A fl „ 



r) 



i i i 



1,1 



4A f 



2' 2' 2' x ' x > (l+A Hl +A J?2 )2' (l + A Hl +A K2 ) 2 



(D.12) 

where by F 2 we denote one of the hypergeometric series in two variables The 
parameters A# a ,a = 1,2 play the role of "dressed" coupling constants and are 
defined through 



A„ 



A 



R a 



Cooe^f 1 + (\ Rl + x R2 y 

From these formulae it is possible to investigate the existence of critical behaviour 
in the same way as we did in section 4 for the one dimensional case. Nevertheless, 
this analysis can be carried out directly from eq. ( p. 10 ). To this respect, we impose 
the condition 



a 



2 D 



;d.i4) 



oo o=l 



For large values of | n | the integral ( p.10 ) is dominated by the small momenta 



k | < e, where e ~ is a cutoff; namely 



tr 



nf^ ink (E^) 1/2 



(D.15) 



and assuming isotropy, i.e A a = A, Va, we obtain the scaling behaviour 

1 



11 



\D+1 



(D.16) 



Away from the condition ( p. 14 ) we shall define a "mass parameter" fi 2 as 



^2 = ^ooe^ _ 2D 

A 



(D.17) 



Because the isotropy, the relation between /x and the sole "dressed" coupling 
constant \ R existing in this case, comes from ( p,13|) and is given by 



DX 



R 



R 



(D.18) 



In terms of /i, the behaviour of the t n for large values of | n | is dictated by 

/ f[ ^ e ink (/i 2 + | k | 2 ) 1/2 (D.19) 



Tfp 
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Moreover, a saddle point evaluation of this integral for fi 2 > yields the expo- 
nential decaying law 



| n 1 2 +i 

Hence, defines the correlation length. 

The constraint to = — 1 gives us an expresion where the density of energy 
and the parameter \i are tied together, namely 



e^ p \ 1/2 r^dh r^dk 



XI Jo 2tt Jo 2n 



/i 2 + 4£sin 2 (^) 

a=l A 



1/2 

(D.2i; 



To investigate the singular behaviour of e^f as a function of /1, it is sufficient to 
consider the small values of the "momenta" k a as we did previously for the t n . In 
this case the integral (p.21|) becomes proportional to the function F D (fi,e) defined 



as 



F D (u, e) = P dk k D ~ l (/i 2 + £; 2 ) 1/2 (D.22) 
Jo 



where e is a cutoff. Studying the singular behaviour of F D (fi,e) in the vecinity of 
the point /i = we get: 

( c D n D+l log(f ) + F r D e9 {^ e) for D odd 
F D (ji,e) = l (D.23) 
[ F D e9 ([i,e) for D even 

where F D eg (fi, e) denotes a regular function at /x = 0. Eq. (|D.23| ) means that only in 



odd dimensions (D — 3, 5, ■ • •) we obtain a singular logarithmic behaviour; that is 



somehow the generalization of the one dimensional result ( 4.38 ). The only difference 
with this latter is the power n D+1 in front of the logarithm, which is expected from 
dimensional grounds. 

Curiously enough, we do not observe singularities in the ground state energy if 
the dimension is even. This fact does not exclude however the existence of discon- 
tinuities in some derivative of the energy as we cross the critical point, just as it 
occurs in a mean field calculation. 

In any case, it is clear that the ansatz (|D . 1|) does not lead for D > 1 to the 
standard ITF model. That supposes, according to the general ideas expressed in 
the introduction of the paper, that we have to return to perturbation theory and 
figure out the terms necessary to improve the ansatz. In D = 2 for example, we 
have neglected in the ansatz the terms involving four a z matrices. We claim here 
that these terms are important if one wants to get closer to the ITF model. As a 
result, we propose the following ansatz for D > 2 



V>({a ri , r2 , r3 })> = exp( £ -f^- £ a^ +ri a? +r2 af +r3 ) | 0) (D.24) 

r'sGA' Zn ri,r 2 ,r 3 j e A 



where n rijr2ir3 is defined by 
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_ J 2 if r Q = v b ^ r c , , 

n ri ,r 2 ,r 3 - | 4 Qtherwise \U-*>) 

Equation ( |D.24j ) includes the ansatz ( |D.1| ) by allowing two r's to coincide. The 
quantities (|D.25| ) accounts for the overcounting in the sum over j in the exponent of 



indices and must satisfy the reflection symmetry property 



( p.24| ). Likewise the parameters a ri)r2)r3 has to be totally symmetric in their three 

(D.26) 



which is the analogue of condition a r = a_ r in the one dimensional case. 
The FP-Hamiltonian associated to (|D.24j ) is 



H FP = ~ XX - eX P(" E ^r^Vi+rrf+rrf+J} (D.27) 
j6A r'seA' 



and can be expanded as 



# FP = -£ FP -E^-E 7^^E^^^^ + (6^b) + ... (D.28) 



jeA r'seA' l n ri,r 2 ,r 3 y- j eA 



Compairing (|D.27| ) and (|D.28j ) and using the same trick of section 4, according 
to which one could replace consistently the cr? by a- } = ±1, we obtain the following 
equations for E FP and A riir2jr3 as functions of a riir2ir3 

EFP = ~^E eX P(~ E "ri.ra.rs^Offri^ra^rs) (D.29) 
cr's r'seA' 

A Sl ,s 2 ,s 3 = -^7E(E°"j (T j+si Cr j+S2 (T j+S3) eX P(- E «r 1 ,r 2 ,r :i O-0O"r 1 O"+r 2 Crr 3 ) (D.30) 
cr's jeA r'seA' 

Naturally the problem which arises now is to perform the sum over the cr's. 
Observe that eq. ( |D.29| ) is a kind of classical partition function of a lattice involving 
couplings between 2 and 4 spins; in some sense, it recalls an 8-vertex model |]J so 
that there are chances to perform the sum, at least in some particular cases. 

Finally, there is the question of inverting eq.( p.30|) in order to express the a's as 



functions of the A's. Notice that the A's could be thought of as a kind of correlators 
among cr's for a partition function as the one of eq.( |D.29| ). 



In summary, the problem of constructing the ground state for the 2D quantum 
ITF model is essentially reduced, by means of the ansatz ([D.24|) , to that of computing 



a 2D nonlocal partition function followed by the inversion of the relation between 
certain correlators and coupling constants. The difficulties of these two tasks are 
considerable but we expect that the FP procedure can throw some light on the 
original unsolved problem. 
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Table captions 



Tablel: Values of the coupling constants and interaction terms in the case of a ITF 
model defined on a finite lattice with N = 4 sites. 
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Figure captions 



Figure la).- The shadow region simbolizes the possible non vanishing entries of 
the matrix p n I . 

Figure lb).- Information obtained up to order v in perturbation theory. 
Figure lc).- Information obtained up to order v using a truncated cluster expan- 
sion. 

Figure Id).- Here we choose one or more representants for each order in pertur- 
bation theory and construct a Hamiltonian which belongs to the same universality 
class of the original Hamiltonian. 

Figure 2).- Comparison of the values for minus the ID energy density in the cases 
exact tising,D=\ (strong solid line), in the first order FP-aproximation ef s f„^ D=1 (solid 

line) and in the second order FP-aproximation ef^^ D=1 (grey line). 
Figure 3).- Action of the operator A on a generic state | I a ) —\ jii, ■ ■ ■ , ji N ). 
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Table 1 



Labels 


Interaction 


Coupling constant 
in the ansatz 


Coupling constant 
in the hamiltonian H FP 


h = (0,0,0,0) 











h = (1,1, o,o) 

h = (0,1, 1,0) 
h = (0,0, 1,1) 
h = (1,0, 0,1) 






A 


/ 6 = (1,0, 1,0) 
/ 7 = (0,1, 0,1) 




a 2 





h = (1,1, 1,1) 




«3 
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